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Abstract: Applying Artmann's formula to the TE-polarized incident 
waves, we theoretically show that the Goos-Hanchen (GH) shifts near the 
angle of the pseudo-Brewster dip of the reflection from a slab of chiral 
metamaterial can be greatly enhanced. The GH shifts are observed for both 
parallel and perpendicular components of the reflected field. In addition, 
it is found that the GH shifts depend not only on the slab thickness and 
the incident angle, but also on the constitutive parameters of the chiral 
medium. In particular, when the incident angle is close to the critical angle 
of total reflection for LCP wave, significant enhancement of the GH shifts 
can be obtained. Finally, the validity of the stationary-phase analysis is 
demonstrated by numerical simulations of a Gaussian-shaped beam. 
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1. Introduction 

The Goos-Hanchen (GH) effect [1, 2], which usually refers to the lateral shift of the totally 
internal reflected from the position predicted by geometrical optics, has been widely analyzed 
both theoretically [3, 4, 5, 6] and experimentally [7, 8, 9]. This phenomenon has already been 
extended to many fields such as acoustics, surface optics, nonlinear optics, and quantum me- 
chanics [10, 11]. Furthermore, with the development of near-field scanning optical microscopy 
and lithography [12], the importance of the GH shift has been paid more and more attention. For 
instance, the GH shifts were found to be large positive or negative for both reflected and trans- 
mitted beams in different media or structures, such as dielectric slabs [13, 14, 15, 16], metal 
surfaces [17, 18], dielectric-chiral surface [19, 20], absorptive media [21, 22], multilayered 
structures [23] and photonic crystals [24]. Recently, the GH shift associated with left-handed 
material has been extensively studied [25, 26, 27, 28, 29, 30] owing to its very unusual prop- 
erties. On the other hand, since Pendry proposed the Swiss roll structure to achieve a chiral 
medium with negative refraction [31], negative chiral medium has been receiving great interest 
[32, 33, 34, 35]. In [36], we have discussed the GH shift of electromagnetic waves incident 
from an ordinary medium into negative chiral medium half space. But, the thickness of chiral 
metamaterial is always finite, and the properties of the GH shift and the reflection or transmis- 
sion will strongly depend on the thickness. Therefore, it is worth discussing the effect of the 
thickness of the chiral metamaterial on the GH shift and the contribution of the thickness and 
the incident angle to the GH shift together This is the main purpose of the present paper 

In this paper, we study the GH shifts of the reflected waves from a chiral metamaterial slab. 
We predict that the GH shift near the angle of the pseudo-Brewster dip from such a slab can 
be large, and both positive and negative lateral shifts are possible. It is also shown that the GH 
shift depends on the thickness of the slab, the incident angle of the wave and the constitutive 
parameters of the chiral medium. Finally, numerical simulations for a Gaussian-shaped beam 
have been performed to confirm the validity of theoretical results. Here, only perpendicularly 
(TE) polarized incident wave is discussed below, and the results for parallel (TM) polarized 
incident wave can be easily obtained in the same way. 

2. Formulation 

The configuration for the chiral slab in a vacuum is shown in Fig. 1. We assume that a linearly 
polarized wave propagating in a homogeneous isotropic dielectric medium is incident at an 
angle 0, upon the surface of a chiral slab with the thickness d. The constitutive relations of the 
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Fig. 1. Schematic diagram of a light beam propagating through the chiral slab placed in 
free space. 

chiral slab are defined by [37] 

D = ££oE - iK^EoHoa, (1) 
B = ^^oH + /Kr^eoMoE, (2) 

where K is the chirality parameter (assumed to be positive in this paper), e and jJ. are the rela- 
tive permittivity and permeability of the chiral medium, respectively (Eq and IM) are the permit- 
tivity and permeability in vacuum). For simplicity, a monochromatic time-harmonic variation 
exp(/a)r) is assumed throughout this paper, but omitted. It is seen that two interfaces of the 
chiral slab are located at z = and z^d. 

The incident electric and magnetic fields of an incident TE wave can be written as 



Ei = £'oi'exp[;^,(— cos 0,z + sin 



(3) 



H; = . / — £■()(— COS 0,i — sin 0,£)exp[;A:,(— COS + sin 
V Mo 



(4) 



where k, = ko = co^Eq^q is the wave number in the air. In a chiral material, an electric or 
magnetic excitation will produce both the electric and magnetic polarizations simultaneously 
[38, 39]. Therefore the reflected wave must be a combination of the perpendicular and parallel 
components in order to satisfy the boundary conditions. According to the propagation direction 
of the reflected wave, the electric and magnetic fields are expressed as 

Er — E()[R iiy + R2i sin OrZ — R21 COS 9rx] e\p[ikr{cos 0,-1 + sin Orx)], (5) 

Hr = A —EQ[R2iy + Riicos 9,-x — Riisin9rz]exp[ik,-{cos 9rZ + sin9,-x)], (6) 
V Mo 

where 9r is the reflected angle, k,- = ki is the wave number of the reflected wave, R 1 1 and R2 1 are 
the coefficients associated with perpendicular and parallel components of the reflected wave, 
respectively. 

Inside the chiral slab, there are two modes of propagation: a right-circularly polarized (RCP) 
wave at phase velocity (o/ki and a left-circularly polarized (LCP) wave at phase velocity to 7^2 ■ 
The wave numbers ki and ^2 have the form [37] 

ku2^ko{^±K). (7) 

Then, the refractive indices of the two eigen-waves in the chiral medium are thus given as 

ni,2 = VeMiK^- (8) 

Recent studies have shown that k > ^/Eji can occur at least at or near the resonant frequency 
of the permittivity of a chiral medium (called chiral nihility [32]), and then backward wave 
will occur to one of the two circularly polarized eigen-waves, making negative refraction in 
the chiral medium possible. When K > ye/i, the refraction index «i = ^/eji+ K will still 
be positive, but the refraction index «2 = \/£M ^ ^ will become negative. Correspondingly, 
negative refraction will arise for LCP wave. Based on the discussion above, in the chiral slab, it 
is assumed that there exist four total waves, two propagating toward the interface z^d and the 
other two propagating toward the interface z = (see Fig. 1). The electric and magnetic fields 
of these waves propagating inside the chiral medium toward the interface z = d axe. written as 

E+ ~Eq[Ay cos 9\x+Aisin9iz + iMy] exp[i^i (— cos 9yz + sin 9\x)\ 

+E(j [A2 cos 92X+A2 sin 92Z — iA2y] exp[;A:2 (— cos 92Z + sin 62-^)] , (9) 
= ;T7"*£'o[^icos 9ix+Ai sin9iz + iAiy]exp[iki{—cos 9iz + sin9ix)] 
—irj^^Eo [A2 cos 92X+A2 sin 92Z — iA2y] exp[/A:2 (— cos 02Z + sin 62^1^)] • (10) 

The total electromagnetic fields of the other two waves propagating inside the chiral medium 
toward the interface z = are expressed as 

= Eq[—Bi cos 9ix + Bi sin + /-Bi3']exp[;^i (cos 9iz + sin9ix)] 
+Eo [—B2 cos 92X + B2 sin 92Z — iB2y\ exp[;A:2 (cos 622 + sin 92x)] , (11) 
= /77^'£'o[— ^1 cos 9ix + B\ sin9\z + iB\y\ exp[/A:i (cos 9\z + sin0i.x:)] 
—ir\ ^^Eq [—B2 cos 92X + B2 sin 92Z — iB2y] exp[/A:2 (cos 92Z + sin 92x)] , (12) 

where rj = -^//i/e is the wave impedance of the chiral medium, A 12 and Bi 2 are the transmitted 
coefficients, and 0i 2 denote the refracted angles of the two eigen-waves in the chiral slab, 
respectively. 



Outside the slab (z > d), the total transmitted wave can be expressed as 



Et = £'o[7"iii'+ 721 sin0,£+ 721 cos 0,1] exp[/^,(— cos Qiz + sin Qtx) 



-EolTwy— T21 sin0,£ — T21 cos 9tx\ exp[;^,(— cos OtZ + sinOtx)], 



(13) 
(14) 



where kt = kj, 0, is the transmitted angle, Ti\ and T21 are coefficients associated with perpen- 
dicular and parallel components of the transmitted wave in the air, respectively. 

The coefficients Rn, R2i,A\2^Bi2,Tii and T21 can be determined by matching the boundary 
conditions at two interfaces z~Q and z~d, and the following matrix can be obtained: 
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(15) 



where 



mil = 

^21 = 



^22 = 



— /(tj cos 0,- + COS 01 ) — cos 0/ — 77 cos 0i 
/(77COS0, — cos0i) cos 0,- — 77 cos 01 

/(tj cos 0,- + cos 02) — COS 0/ — 77 COS 02 
— /(?] COS 0; — COS 02) COS 0,- — 77 COS 02 

/(-77 COS 0/ + COS 01 )e-'('^-c-'*^u)^ (- cos 0/ + 77 cos 0i )e-'(^';-'^ic)'' 
/(77 COS 0,- + COS 01 )e-'^'''~+'''z)d (77 cos 01 + COS 0,)e-'(<:'c+'^l:)'' 

/(77 COS 0; - COS 02)e-'('^«-'^2z)<' (_ cOS 0; + 77 COS 02)e-'('^':-'^2c)^/ 

-/(77 COS 0,- + COS 02)e-'(*^«+*:2j)<' cos 02 + COS 0,)e-'(*^':+*^2;)'' 



(16) 
(17) 
(18) 
(19) 



The analytic solutions to the above matrix can be obtained after some lengthy mathematic 
manipulations, but the final results are too complicated to reproduce here. Now we would like 
to study the GH shifts of the reflected beams from a slab of negative chiral medium. It is well 
known that the lateral shift of the reflected beam has the definition 



A = d(i>/dkx, 



(20) 



which was proposed by Artmann using the stationary -phase method [3]. <t> is the phase differ- 
ence between the reflected and incident waves. For a chiral medium, the reflected field contains 
both parallel and perpendicular components, these field components can to first order be repre- 
sented as two separate reflected beams, each with its own magnitude and lateral shift. Then the 
GH lateral shifts for perpendicular component ( 7? 1 1 in Eq. (15)) and parallel component {R21 in 
Eq. (15)) of the reflected wave can be calculated by the stationary-phase approach as a function 
of the angle of incidence. 



3. Results and discussion 

In this section, numerical simulations associated with the GH lateral shifts for a negative chiral 
slab will be presented and compared with a conventional chiral slab. 

In our calculation, we consider two types of chiral slabs: (1) a positive chiral slab (e.g., 
£ = 0.64, jj. = I, K = 0.4) whose refraction indices of RCP and LCP waves are both positive, 
i.e, a conventional chiral slab; (2) a negative chiral slab (e.g., JC = 1 .4 and the other parameters 
unchanged) whose refraction index of LCP wave is 112 ~ —0.6. The working frequency here 
is © = 27r X 10 GHz, and the thickness d of the chiral slab is 1.5A. Figs. 2(a) and 2(b) show 
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Fig. 2. The dependence of (a,b) the lateral shifts and (c,d) the phases of the reflection 
coefficients for (a,c) perpendicular and (b,d) parallel components on the incident angle in 
the presence of different chiral slabs. The insets of (a) and (b) show the absolute values of 
perpendicular and parallel reflection coefficients, respectively. Solid line and dashed line 
correspond to positive chiral slab and negative chiral slab, respectively. 

the dependence of the GH shifts of the reflected beam upon the incident angle for both per- 
pendicular and parallel components from the slab of Type- 1 and Type-2. In Figs. 2(a) and 2(b), 
both the insets show, respectively, the absolute values of perpendicular and parallel reflection 
coefficients for two different types. It is easily found that there is a dip in each reflection curve, 
at which \R\ reaches the minimum (close to zero). The coiTesponding incident angle is defined 
as the pseudo-Brewster angle. From Figs. 2(a) and 2(b), we clearly see that the behaviors of the 
GH shifts for perpendicular and parallel components are similar, and the shift will be greatly 
enhanced near the angle of pseudo-Brewster dip. For the negative chiral slab, the lateral shifts 
of both reflected components might be large negative near the dip, and be small positive values 
under other incident angles. This phenomenon can be easily explained in terms of the change 
of phase. The phases of reflected parallel and perpendicular components as a function of the 
incident angle are plotted in Figs. 2(c) and 2(d) for negative and positive chiral slabs. Near the 
angle of the dip, the phase of reflection experiences a distinct sharp variation, which decreases 



quickly for the chiral negative slab. As a result, Artmann's formula Eq. (16) leads to a large 
negative GH shift and small positive lateral shift. In contrast to chiral negative slab, for the 
positive chiral slab, both components have positive lateral shifts near the angle of dip, and then 
experience small negative shifts over other incident angles. The dependence of GH shifts on 
incident angle for negative chiral slab is opposite to that for positive chiral slab. 
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Fig. 3. (a) The dependence of the lateral shifts for reflected perpendicular and parallel 
components on the incident angle for a typical negative chiral slab, (b) The absolute values 
of perpendicular and parallel reflection coefficients. The inset of (b) is the phase of both 
reflection coefficients. 

Figure 3 shows a typical dependence of the GH shifts of reflected perpendicular and parallel 
components on the incident angle from a negative chiral slab with a large chirality parameter. 
Assume that the chiral slab has a chirality K" = 2.0, while other parameters remain the same 
as before. Figure 3(b) shows the absolute values of the reflection coefficients for perpendicular 
and parallel components as a function of the incident angle. It can be seen that there is only 
one dip in the perpendicular reflection curve, at which the reflection coefficient reaches a min- 
imum magnitude and the phase difference increases monotonically as a function of incident 
angle (see the inset of Fig. 3(b)). Thus the GH shift of the reflected perpendicular component 
has a positive peak near the angle of the dip. In contrast, there are two dips in the parallel 
reflection component, where the absolute values of the reflection coefficient are very close to 
zero. Meanwhile, from the inset of Fig. 3(b), the corresponding phase in the vicinity of these 
two dips monotonically decreases quickly. Therefore, it indicates that the shifts of the parallel 
component can be greatly enhanced to be large negative near certain angles where the phase 
decreases. The shift enhancement can be an order in magnitude greater than the wavelength. 

In what follows, we turn to discuss the lateral shift as a function of the thickness of the slab 
under different incident angles. First, we consider the GH shifts of the reflected beam from a 
negative chiral slab. As an example, the dependence of lateral shifts on the slab thickness at 
different incident angles is presented in Fig. 4. The parameters are taken as follows: e — 0.64, 
jj. = I, K ~ 1.3, implying that two refraction indices corresponding to RCP and LCP waves are 
«i =2.1 and «2 = —0.5. It is obvious that there exists a critical angle at Oc = 30° for the LCP 




Fig. 4. The dependence of lateral shift on the thickness of the chiral negative slab for differ- 
ent incident angles, (a) and (c) Reflected perpendicular component; (b) and (d) Reflected 
parallel component. 

wave, but it does not correspond to the true total reflection. Under such critical situation, only 
LCP wave in the chiral slab becomes evanescent wave when the incident angle exceeds this 
critical angle, while the RCP wave will still propagate through the chiral slab. Correspondingly, 
this incident angle is defined as pseudo-Critical angle. However, the true total internal reflection 
will never arise under these parameters. In such a case, when the incident angle is smaller 
than the critical angle (9c = 30°), it is shown in Fig. 4(a) that the GH shift of the reflected 
perpendicular component could reach large negative or positive values, but there never exists 
any periodic fluctuation of the lateral shift with respect to the thickness. Similar phenomenon 
can be found from the shift of the parallel component (see Fig. 4(b)). Their large negative 
enhancement of the lateral shifts correspond to the dips of the reflection coefficients, which 
follows previous discussions. However, when the incident angle is equal to the critical angle of 
the LCP wave, as shown in Figs. 4(c) and 4(d), the lateral shifts of both reflected components are 
always positive and increasing with respect to the thickness of the slab. This can be explained 
by the fluctuation of the reflection coefficient with the slab thickness. However when the slab is 



much thicker, the strong fluctuation of GH shift can be attained, as shown in the insets of Fig. 
4. And when the incident angle further increases above the pseudo-Critical angle, the small 
fluctuation of the lateral shift against the slab thickness is observed, as shown in Figs. 4(c) and 
4(d). 




Fig. 5. Dependence of the GH shift on the thickness of an invisible medium slab for differ- 
ent incident angles, e = 0.8, fl = 0.8, and K = 0.2. 

Apart from the above-mentioned the GH shift of the negative chiral slab, we also consider 
the lateral shift of the other chiral metamaterial slab. Here, we set the parameters of the chiral 
medium as e = 0.8, ji = 0.8, and K = 0.2. In this case, the wavenumber matching condition ki = 
ko and the wave impedance matching condition 77 = 770 are satisfied simultaneously, the RCP 
wave is transmitted straight through the chiral medium without either reflection or refraction 
at any incident angle. This medium is invisible for RCP wave [42]. While the LCP wave can 
be refracted and reflected, or totally reflected. This unusual phenomenon can be physically 
understood as a destructive interference of electric and magnetic responses, due to the mixing 
through the chirality parameter For ^2 ~ 0.6^o, Snell's equation for LCP wave is expressed as 
sin 0, = 0.6sin 62, so the critical angle for LCP wave is 9c = arcsinO.6 ~ 37°. Therefore, LCP 
wave is totally reflected with the incident angle greater than 37°. This implies that we can divide 
the incident wave into LCP and RCP waves. The property of GH shifts for this case is shown in 
Fig. 5. Fig. 5(a) shows the perpendicular and the parallel reflected GH shifts as a function of the 
slab thickness when the incident angle is equal to the critical angle of the LCP wave (Be = 37°). 
Calculation under these conditions shows that both reflected components have the same GH 
shift. This is because the reflected wave only has LCP wave, and the RCP wave contributes to 
the transmitted wave. Hence the perpendicular and the parallel reflection coefficients have the 
same absolute value, while their phases are different. From Fig. 5, it is straightforward to see 
that the calculated lateral shifts for this case are negative. Moreover, when the incident angle is 
close to the critical angle of the LCP wave, the lateral shifts are large and increase as the slab 
thickness increases. While the incident angle is greater than 9c, we can see that, as the increase 
of the slab thickness, the lateral shifts will increase quickly and then gradually approach to an 
asymptotic negative value. Therefore, the GH shift for the thick slab in the present case depends 



on the incident angle more than the slab thickness. 
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Fig. 6. Dependence of the GH shift on the incident angle. The theoretical result is shown by 
the line; the numerical results are shown by solid scatters (for wq = lOA) and open scatters 
(wfl = lOOA), all the other optical parameters are the same as in Fig. 3(a). 



4. Goos-Hanchen shift of a Gaussian-shaped incident beam 

To show the validity of the above stationary-phase analysis, numerical calculations are per- 
formed, which confirm our theoretical results. In the numerical simulation, an incident 
Gaussian-shaped beam is assumed,£',(x,z = 0) = exp(— x^/2w^ + ikj^Qx), which has the Fourier 
integral of the following form: 



£,(x,z = 0)= J A{k,)expiik,x)dk, (21) 

where Wx = wo sec 0,, wois the beam width at the waist, and the amplitude angular-spectrum 
distribution is Gaussian, Alk-^) — expf— — kx-o)^]. So the electric fields (£1 and ) 

yLTl ^ ' II 

of the reflected beam, determined from the transformation of the incident beam, can be written 
as 



£'^(x,z = 0) = y RnA{kx)&y.p{ik„x)dky:, (22) 
El,(x,z = 0)= [ R2iA{kx)exp{ikxx)dkx. (23) 



IS 

z=0 



The calculated beam shift can be obtained by finding the location where [is^l, g or E^^ 
maximal [14]. 

Calculations show that the stationary-phase approximation for the GH shift is in good agree- 
ment with the numerical result. As an example. Fig. 6 shows the numerical calculation results 
of curves in Fig. 3(a). The incident beam width is chosen to be wo = lOA and wq = lOOA. For 
comparison, both the numerical and theoretical results are shown in Fig. 6. The peaks of the 
numerical shifts for the perpendicular reflected field are about 6.4A for wq ~ lOX and 7.4A 
for Wo = lOOA, and for the parallel reflected field are about — 6.1A (dip I), — 16A (dip II) for 
wo = lOA, and -9.8A (dip I), -29.0A (dip II) for wq = lOOA. The peaks of the theoretical 
shifts are about 7.39A for the perpendicular field and — 10.07A ( — 29.32A) for dip I (dip II) of 
the parallel field. It is noted that the discrepancy between theoretical and numerical results is 
due to the distortion of the reflected beam, especially when the waist of the incident beam is 
narrow [40, 41]. The further numerical simulation shows that the wider the incident beam is, 
the less the discrepancy is, and the closer to the stationary-phase result the numerical result is. 



5. Conclusion 

To summary, an investigation on the GH shifts of both reflected parallel and perpendicular 
components for chiral metamaterial slab has been done by using the stationary-phase approach. 
It shows that the GH shift of the reflected perpendicular components can be greatly enhanced 
to be large negative as well as positive near the dip of the reflection curve, at which \R\ reaches 
a minimum. Similar behavior will be found for the parallel component. These results obtained 
from the negative chiral slab are opposite to those from the conventional chiral slab. In addition, 
at a given incident angle, the dependence of the GH shift on the slab thickness for negative chiral 
slabs has also been studied. It is shown that, when the incident angle is equal to the critical 
angle of the LCP wave, the GH shifts of both reflected components oscillate with respect to the 
thickness of the slab, and its overall tendency is increasing along with the small slab thickness, 
while for a much thicker slab, the fluctuation of the lateral shift with the slab thickness becomes 
more obvious. However, greatly enhanced GH shifts for a smaller thickness slab can be obtained 
when the incident angle is smaller than the critical angle. Meanwhile, we calculated the GH 
shift of an invisible chiral medium for RCP wave. Finally, in order to demonstrate the validity 
of the stationary -phase approach, numerical simulations are made for a Gaussian-shaped beam. 
Calculation results show that the wider the waist of the incident beam is, the closer to the 
stationary-phase result it is. 
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